Nonequilibrium kinetics of a disordered Luttinger liquid 
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We develop a kinetic theory for strongly correlated disordered one- dimensional electron systems 
out of equilibrium, within the Luttinger liquid model. In the absence of inhomogeneities, the model 
exhibits no relaxation to equilibrium. We derive kinetic equations for electron and plasmon distribu- 
tion functions in the presence of impurities and calculate the equilibration rate 7_b. Remarkably, for 
not too low temperature and bias voltage, is given by the elastic backscattering rate, independent 
of the strength of electron-electron interaction, temperature, and bias. 
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Introduction. — Interacting electrons in one dimension 
(ID) |l[ have become a focus of interest in nanophysics 
from both the fundamental and applied perspectives. Re- 
cent technological advances have made it possible to sys- 
tematically study the transport properties of a variety 
of ultranarrow wires; in particular, carbon nanotubes [2J 
and semiconductor nanowires [3| . From the fundamental 
point of view, much of the fascination with physics of ID 
systems is driven by the fact that electron-electron (e-e) 
interactions in ID geometry are qualitatively significant, 
transforming the electron gas into a Luttinger liquid (LL) 
[l[ . This strongly correlated state of matter is commonly 
described in terms of bosonic elementary excitations. 

A conceptually nontrivial aspect of the non-Fermi liq- 
uid nature of a LL concerns its behavior at nonequilib- 
rium, e.g., when a finite bias voltage is applied to the 
wire (for a recent experiment see Ref. Q). A homoge- 
neous LL is completely integrable and as such does not 
exhibit any relaxation to equilibrium: an excited state 
will never decay to the state characterized by tempera- 
ture T. Quite remarkably, a finite quantum lifetime of 
fermionic excitations due to e-e interactions in a homo- 
geneous LL does not translate into any inelastic e-e scat- 
tering [j|: the allowed energy transfer is exactly zero. 
This is in stark contrast to electron liquids in higher di- 
mensions, where the characteristic energy transfer is T. 
Relaxation to equilibrium due to e-e collisions in a LL 
is thus only possible if momentum conservation is bro- 
ken by inhomogeneities. Of central importance is there- 
fore the question — essentially unanswered — of how the 
equilibration in a LL occurs in the presence of a random 
backscattering potential. This is the subject of this work. 

So far, advances in dealing with a LL off equilibrium 
have been focused on a "mechanical" approach, i.e., on 
solving the equations of motion "as exactly as possible" . 
Efforts, based on the bosonization approach, have been 
centered around the nonlinear conductance of a LL con- 
taining a single compact scatterer Q. However, the line 
of research relying on the exact integrability cannot pos- 
sibly be much extended beyond the single-scatterer case 
(for two tunneling barriers in a LL, a nonequilibrium dis- 



tribution of plasmons was studied by means of the mas- 
ter equation, see Ref. [7[). Also, importantly, the conven- 
tional bosonization [l|] is designed for equilibrium bound- 
ary conditions. An alternative is to preserve dynamics of 
both bosonic and fermionic degrees of freedom [1, @] , as 
it is done in the functional bosonization approach (see 
Ref. for a review). 

Our purpose here is to develop a kinetic approach 
to nonequilibrium phenomena in a disordered LL, by 
formulating kinetic equations for distribution junctions. 
Within this approach, one has to introduce the distri- 
bution functions of not only bosonic but also fermionic 
excitations, similarly to higher dimensions ll|, LL2|. Our 



main result is a set of kinetic equations which describe 
(i) inelastic e-e scattering, mediated by virtual plasmons, 
and (ii) creation/annihilation of real plasmons; both pro- 
cesses being only triggered by scattering off disorder. We 
calculate a key quantity in nonequilibrium problems: the 
energy relaxation rate je- In a remarkable departure 
from Fermi liquids, in a LL at not too low T turns 
out to be given by the elastic scattering rate 7. 

Effective action. — We study interacting electrons in a 
single-channel disordered quantum wire within the LL 
model [lj : the electron dispersion relation is linear (with 
the velocity Dp) and interactions yield only forward e-e 
scattering (characterized by the dimensionless constant 
a = V{/2ttvf, where V{ is the zero-momentum Fourier 
component of the interaction potential). We consider 
both spinless (rj = 1) and spinful (77 = 2) models. Our ap- 
proach is based on the "quasiclassical" real-time electron 
Green's function at coinciding spatial points g(x, tx^t^) 
[HI, widely used in the nonequilibrium theory of meso- 
scopic transport. The "hat" means that g is a matrix 
in the Keldysh, chirality (/1 = ± for right/left- moving 
electrons) and (possibly) spin (s =f,J,) spaces. In what 
follows we use the Pauli matrices r, cr, and s that act in 
the chirality, Keldysh, and spin spaces, respectively. 

The quasiclassical Green's function satisfies the con- 
dition g o g = 5(ti — £2); where the "dot" denotes the 
convolution in the full (chirality x spin x Keldysh x time) 
space. This constraint enables us to describe the wire 



by the action that reproduces the averaged over disor- 
der equation of motion for g as its saddle point. The 
way to derive such an action is similar to that in the 
ballistic sigma-model in higher dimensions 14(. To ac- 



count for the Coulomb interaction, one introduces the 
Hubbard-Stratonovich field <fi(x,t) = 4>i + Gxfyi on the 
Keldysh contour [lij], where the "classical" and "quan- 
tum" fields 0i,2 are diagonal matricies in the chirality 
and spin spaces. Then the action takes the form Q 
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Here 70 is the bare elastic rate of backscattering off a ran- 
dom static potential, go = diag (g^ , —g^) corresponds to 
the saddle point of the action of the noninteracting prob- 
lem, and g — Tc/qT^ 1 = diag (g + , — g~~). The unitary 
transformation T (diagonal in the chirality x spin space) 
parametrizes fluctuations around go due to fluctuations 
of <j)(x, t). The interaction of fermions of the same chiral- 
ity and spin is included in the shift of the Fermi velocity 
v* F = v F + Vf/2n [5]. The last term in Eq. p with 

actions. A similar (replicated imaginary-time) action for 
a disordered LL at equilibrium was proposed in Ref . [l6| . 

Quantum kinetic equations. — Starting from the effec- 
tive Keldysh action (fl]), we use the standard procedure 
[l3L H3] to derive the kinetic equations. We proceed at 
one-loop order with respect to the effective interaction, 
which is equivalent to the "dirty random-phase approx- 
imation" (dRPA) [5[. The one-loop derivation is con- 
trolled by the parameters 7/ max{T, ell} <C 1 and a< 1, 
which is assumed in the rest of the paper (U is the bias 
voltage). We also disregard the localization effects [HI. [lij]. 
For a wire of length L > Vp/j, this limits the applicabil- 
ity of what follows to max{T, ell} > 7\ = "//a 3 ~ v . 

Within the dRPA, we expand the action (TT]) in fast 
(on a scale of the relaxation time for the kinetic equa- 
tion) quadratic fluctuations of (f> around the "slow" 
semiclassical saddle point g^ ow = {<J Z + cr + )5(ti — 
£2) — 2a+f fJ, (x,ti,t2). The electron distribution func- 
tion fV"(x,t) at given energy e is defined via the Wigner- 
transform of / M (x, tx, £2) and is yet to be found from 
the kinetic equations. In this way we obtain the Gaus- 
sian action with the propagator V = (V^ 1 — ft) -1 . 
Here n = o~ x {dtD a x — X)j 2txv* f is the polarization op- 
erator, with D being the electron-hole (e-h) propagator 
damped by disorder. The retarded part of D is given by 
D^ 1 = —i(u) — T z v* F q) + 7(1 — t x )/2, while the kinetic 
part Dk is expressed (see below) in terms of fj? via 



FIG. 1: Scattering processes corresponding to the collision 
kernels fC(ui) and C(ui). Thin wavy lines: bare interaction 
Vf. Solid wavy lines: dRPA interaction propagator V>. The 
triangles (a) and rhombi (♦) mean the Keldysh part of the 
electron and boson Green's functions, respectively. 

out the fast fluctuations in the equation of motion for 
g(x,ti,t2) with the dRPA action, which yields [l!^] 



7 + 7inol 



(/f-/r")+st t ( £ ). (3) 



The collision integral St^_ & describes inelastic electron 
scattering due to interaction with the bosonic bath, 



m - (4) 

while 7inoi accounts for the additional [191 ] renormaliza- 
tion of the static disorder due to the inelastic scattering, 



7inol 



= - / <kjl+-(u)(l- / e - w +/e +w ) 



(5) 



Here f e = (/+ + f~)/2, and I^{uj) is the rate of emis- 
sion of energy per unit interval of u>, accompanied by 
scattering /i — > i>, which is given by 



= ! / ^V^(to,q)ReD%(to,q) 

TT J ZTT 
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To derive the kinetic equation for electrons, we average 



In the integrand ([6]), there are four poles, q ~ ±cj(1 ± 
i r y/2u) /vp, inherited from Re Dr, which are only slightly 
damped by disorder in the limit w»7. They correspond 
to e-h pair excitations described by the renormalized 
Fermi velocity v F . Four more poles q ~ ±uj{\±i / y/2w)/u, 
associated with the "greater" part of the effective inter- 
action with parallel spins V> 11, correspond to the collec- 
tive plasmon mode of the clean LL, moving with velocity 
u = vp(l + r\a) x l 2 . In the spinful model, V contains 
an extra mode — spinon — propagating with velocity i>f 
[Lsl ]. Importantly, the e-h and collective excitations at 
uj ^S> T\ are well resolved from each other and should be 
treated separately, while in the opposite limit, w C Ti, 
the disorder-induced quantum uncertainty makes them 
indistinguishable. We now proceed by relating V^t to 
the distribution functions. 

Small energy transfer, u) <C T±. — In this limit, the 
Keldysh, retarded, and advanced parts of interaction sat- 
isfy Vk — VrIIkVa- Then the general result for the 
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emission rate can be represented in the form I^ v {uj) = 
J] Q/ 3 w/C^(w)7V" /3 , where the collision kernel /C^(w) 
describes inelastic spin-conserving electron scattering 
fj, — > v with energy transfer oj to the electron and hole 
having the chirality a and (3. For the total collision kernel 
/C = | £ M „ ICfZ we obtain 



^ [IV 



[IV 



(7) 



In the spinless model, the kernel fCiui) is determined by 
processes described by the Feynman diagrams in Fig. QJ, 
at 7 < w < aTi and Fig. [TJd at aTi < w < T x , In the 
spinful case, collisions between electrons of the same chi- 
rality but opposite spin (Fig. [I];), give the main contribu- 
tion to K,(uj) for all w < Tj. The frequency dependence 
and the asymptotes of JC(u>) are shown in Fig. [2] 

Large energy transfer, oj 3> T\. — In this case, we con- 
sider the contributions to l» v = I£ v + 1% V + 1^ from the 
collective (plasmon/spinon) and e-h poles separately, see 
Figs. [T}I-[Tf. Because of the splitting of the e-h, plasmon, 
and spinon poles, the collision integral would be nonlo- 
cal [llj if expressed solely in terms of the e-h distribution 
functions N£ v [see Eq. ((2)]. By introducing the plas- 
mon/spinon distribution functions n b (u>) we can express 
the emission rate of the bosonic excitations (b = p, s) in 
the local form (v p = u, v s = vf): 



1 b 



(8) 



where the collision kernels C have the scaling form C a b — 
("f/2(jj 2 ){v F /vb) A£ v c and the ^4-factors read 

A^ >p = [ V (v F -v F )+v F ±u]/r 1 v Fl (9) 



A 



-(v* F - v F )/v F , A*£ s = (r)- 1)6^6^. 



In turn, n b (ui) satisfies the kinetic equation 



a/3 

(10) 

where 7;, = "f{vp/v F ) and we used the relation D^f ~ 
2Re£>^"(l + 2Ng v ). This kinetic equation describes 
the decay/creation of plasmons and spinons in/from e- 
h pairs. Finally, the e-h pole in the emission rate ^ 
gives = 8 liV {')/w)N^ ,i . Note that the collision 

kernels C\ b ~ C e h = 7/^ 2 do not contain a as a small 
parameter. These most efficient relaxation processes in- 
volve electrons and bosons within the same chiral branch: 
their rates are resonantly enhanced since v b and v* F are 
close to each other. As a result, a appears in the combi- 



1. It can also be shown that the 



nation a 2 /\v* F — v b \ 

total kernel, L = £ eh + (1/2) X^ Q)i) ££" 6 , is equal to the 
asymptotic value of K,(uj) [Eq. Qj] in the limit uj >■ T\. 

Energy relaxation rate. — We now turn to the energy- 
relaxation rate je in the limit of weak nonequilibrium. 



J3 
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FIG. 2: Frequency dependence of the collision kernels for the 
spinless and spinful models [Eq. 0] for a = 0.1. 



f„ 0.5 




FIG. 3: Electron distribution function in the middle of a 
quantum wire of length L biased by ell = 40T, for T = Ti/2. 
The solid lines correspond to "/td = 0.25, 0.5, 1, 2, 4, and 
8, where td = L/vp and 7 is the elastic backscattering rate. 
The dashed curve shows the limiting Fermi distribution. 

Linearizing the electron and boson kinetic equations, we 
estimate 7^ from the total [20( collision kernel (Fig. [2]) 



IE' 



duu 2 K,{uj) ~ T 2 K,{T). 



(11) 



The characteristic ui for "/p is of order T. Since K,(lj) ~ 
27/w 2 for w Ti, 7e at T >■ T\ does not depend on a 
and is given by the backscattering rate: 



7b ~ 7 



(12) 



(a enters only through the condition on T, so that the 
result is also valid for the long-range Coulomb interac- 
tion). At T < Ti, Eq. (fl"2j) is only valid for short wires 
of length L < uf/7 — otherwise the system is localized — 
and shows that a full equilibration then has no time to 
develop. Note that jp for spinful and spinless electrons 
turn out to be parametrically the same, in contrast to 
the weak- localization phase-relaxation rate 7^ 18]. 

In the limit of strong nonequilibrium, we solve Eqs. (|3j). 
(fl"0|) numerically to obtain the distribution function of 
electrons / c in a wire biased by a voltage U ^> T/e, 
where T is the temperature in the leads, as a function of 
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the distance to the contacts (/ e shows up directly in tun- 
neling spectroscopy [41; fo r experiments on multi-channel 
wires, see, e.g., Ref. [2J). The result, shown in Fig. El 
confirms the estimate (|T2"j) : the scale of Try = L/v* F on 
which / e equilibrates is 7 _1 , despite a -C 1. For small a, 
the inelastic processes involving opposite chiralites (£~t + 
and C\_~), as well as the backscattering of plasmons on 



the boundaries |22l, |23j, |24j, [25j, can be neglected: the 
curves in Fig. [3] are thus a- independent. At rrj ^> 7 _1 , 
/ e approaches the Fermi distribution with the tempera- 
ture T e = V3eU/4n. 

Discussion. — As follows from Eqs. (TT2"|) and (TIT)]) , the 
thcrmalization of electrons occurs on the same time scale 
as the lifetime of bosons. This elucidates a conceptually 
important point: boson decay is a source of the inelas- 
tic relaxation of electrons. Indeed, in the homogeneous 
case, the combination ReDImP in Eq. ([7]) yields zero 
energy transfer in the e-e scattering [5j. However, any 
plasmon scattering broadens the peak in ImV", thus al- 
lowing for a finite transfer (even if D remains free — no 
electron backscattering) [2J]. This is true, in particular, 
in an inhomogeneous LL without impurities but with a 
nonuniform interaction a (a;) [22L I23L EBjl . 

In a disordered LL, the spectral function of dRPA 
bosons is characterized by the rate 75 ~ 7 due to 
the elastic scattering off impurities One sees from 
Eq. fTUJ) that out of equilibrium [when n^{uj) ^ N^] 
the boson scattering is represented entirely as the cre- 
ation/annihilation of e-h pairs. As a result, the inverse 
process — the inelastic electron scattering due to the emis- 
sion/absorption of bosonic excitations — is characterized 
by the same rate, je ~ 7- Impurities induce also the 
anharmonic decay of plasmons, as well as their inelas- 
tic scattering on each other. These processes have been 
neglected as higher-loop corrections to the dRPA: their 
rates are much smaller than je- It is the latter that gives 
the thermalization rate for the full bosonic distribution 
function n(u>) in Eq. (|10p due to the coupling (through 
Nu) to the fermionic function f e . 

Conclusion. — We have formulated the analytical 
framework for disordered LLs out of equilibrium, based 
on the kinetic equations for the boson and electron distri- 
bution functions. We have found the equilibration rate, 
which, remarkably, coincides with the elastic scattering 
rate. The kinetic approach developed here is particularly 
convenient for studying heat transport and current noise 
in strongly-correlated disordered ID systems [26]. 
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